Auger decay, Spin-exchange, and their connection to Bose-Einstein condensation of 

excitons in CU2O 
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In view of the recent experiments of O'Hara et al. |Q on excitons in CU2O, we examine the in- 
terconversion between the angular-momentum triplet-state excitons and the angular-momentum 
singlet-state excitons by a spin-exchange process which has been overlooked in the past. We es- 
timate the rate of this particle-conserving mechanism and find a substantially higher value than 
the Auger process considered so far. Based on this idea, we give a possible explanation of the re- 
cent experimental observations, and make certain predictions, with the most important being that 
the singlet-state excitons in CU2O is a very serious candidate for exhibiting the phenomenon of 
Bose-Einstein condensation. 
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Bose-Einstein condensation has been the subject of 
numerous theoretical and experimental studies. While in 
the recent years many of these studies have focused on 
the Bose-Einstein condensation of trapped alkah atoms 
1^, another possible candidate which can undergo this 
second-order phase transition is the exciton gas in semi- 
conducting materials 

Excitons are much like positronium atoms. They are 
bound states which form between electrons and holes in 
a semiconductor, after the electrons get excited from the 
conduction band to the valence band, usually by some 
laser field. Since excitons consist of two fermions, in the 
limit where their separation is much larger than their 
Bohr radius, they are expected to behave like bosons. 
Many experiments have been performed with excitons in 
CU2O because of the many advantages of this material: 
It has a direct, but dipole- forbidden gap, which makes 
the lifetime of excitons rather long, it has isotropic ef- 
fective electron and hole masses, it does not form bound 
states, biexcitons, or an electron-hole liquid, and finally 
the exciton binding energy is quite large. 

The traditional way of observing the kinetic energy 
distribution of excitons is to look at the recombination 
spectrum, and more specifically the optical-phonon as- 
sisted lines. Since the optical phonons have a very weak 
dispersion relation and since the transition matrix ele- 
ment does not depend on the exciton momentum, the en- 
ergy distribution of the emitted photons essentially gives 
the kinetic energy distribution of the excitons. Many 
experiments [^,^ have demonstrated that excitons do in- 
deed obey Bose-Einstein statistics in the limit of high 
enough densities and low enough temperatures, with the 
luminescence spectrum fitting very accurately to Bose- 
Einstein distributions. This fitting procedure gives the 
temperature and the chemical potential of the gas, since 
these are essentially independent parameters. A crucial 
assumption underlying this procedure, is that very fre- 
quent collisions between the excitons bring the gas to 
a quasi equilibrium, with some time-dependent chemical 



potential and temperature, which in general differs from 
the lattice temperature that is kept very low, below 5 K. 
The typical effective temperature of the exciton gas is on 
the order of 10 up to 100 K. Knowing the temperature 
and the chemical potential, one can deduce the particle 
density, assuming an ideal Bose gas with an experimen- 
tally known total exciton mass. The densities turn out 
to be on the order of 10"'^* cm^"^ from this method. Fol- 
lowing this approach Snoke et al. |^,^ observed that the 
triplet-state (ortho)excitons do not Bose condense, but 
move along lines parallel and closely to the critical one, 
which are adiabats, i.e., along lines with constant entropy 
per particle. This effect has been examined theoretically 
in Ref. ||^ , and has been attributed to a competition be- 
tween the acoustic-phonon cooling of the exciton gas, and 
an Auger heating mechanism which prevents the Bose- 
Einstein condensation of orthoexcitons. Lin and Wolfe 
have also reported in Ref. |7| this tendency of orthoexci- 
tons to move along adiabats, but, most importantly, have 
observed evidence for the Bose-Einstein condensation of 
the angular- momentum singlet state (para)excitons 

However, recently O'Hara et al. have developed an- 
other method for estimating the density of the excitons 
By calibrating their photon detector, they have eval- 
uated the number of photons that are being emitted, thus 
determining the number of orthoexcitons inside the crys- 
tal. They have also estimated the volume of the exciton 
gas by knowing the surface of the area of the laser light 
that creates the excitons, and with the assumption that 
the exciton gas has a typical depth inside the crystal, 
which is on the order of the absorption length of the laser 
light. Dividing the exciton number by the volume, they 
have found that the average density of the orthoexciton 



gas is two orders of magnitude smaller 
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than the one they estimate by the spectroscopic method. 
This implies that the gas should be completely classical, 
without showing any kind of quantum degeneracy. 

Therefore one is confronted with a paradox, since from 
the one point of view the spectra can be fitted very ac- 
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curately to Bose-Einstein distributions, but on the other 
hand the densities seem to be quite lower than those one 
gets by this fit, and certainly much lower than the region 
where one would observe Bose statistics. Furthermore, 
Refs. [p], p^|jTl| ] suggest that such low exciton densities 
might be due to a very effective Auger decay mechanism, 
where two excitons collide, the one recombines, trans- 
ferring its energy to the other, which ionizes. The im- 
plication is that this Auger mechanism, which does not 
conserve the total number of excitons prevents the onset 
of Bose-Einstein condensation. Up to now the Auger pro- 
cess was thought to provide the only relatively fast chan- 
nel for excitons to get destroyed However, in view 
of the very long intrinsic radiative lifetime of orthoexci- 
tons reported recently in Refs. jlj, the Auger decay rate 
should have a giant value, exceeding by three orders of 
magnitude the value calculated in Ref. ||Tl[] . 

In this study we examine another mechanism [T^ , in 
which two orthoexcitons with opposite collide, where 
J is the total angular momentum of each exciton, ex- 
changing their electrons or holes in the process, giving 
two paraexcitons in the final state. We make an esti- 
mate of the rate for this spin-exchange process, and find 
that it is rather high. Based on the result of our calcula- 
tion, we propose that this ortho-to-para interconversion 
mechanism is actually the dominant one in the experi- 
mental conditions used so far. Strong experimental ev- 
idence that our argument is true provide Fig. 2 of the 
first paper in Refs. Fig. 5 of the second paper in Refs. 
0, and Fig. 4(a) of Ref. 0, where for late times it is 
clear that there is very slow decay of the paraexciton 
number. Consideration of this process removes the con- 
tradiction described above between the two methods that 
give the exciton density. An important conclusion drawn 
from this new scenario is that the Bose-Einstein conden- 
sation of paraexcitons is probable, since at late times they 
should form a cold and relatively dense gas. The spin- 
exchange mechanism, even if it converts one species into 
the other, conserves the total number of excitons. Since 
the orthoexcitons lie higher in energy than the paraexci- 
tons due to the exchange interaction by an amount Ai?, 
the interconversion process also transfers energy to the 
exciton gas. We also note that our mechanism explains 
the observed sublinear dependence of the orthoexciton 
number jl] that is generated by the laser pulse as func- 
tion of the laser power . On the other hand, under 
extreme pumping conditions in the band-to-band region, 
we have a highly noncquilibrium system at early times 
and one can think that Auger processes between the free 
carriers are also effective. 

Let us start by making an order of magnitude esti- 
mate of the decay rate To.p of the ortho to paraexciton 
conversion process. We consider the process of two or- 
thoexcitons with momenta K and P and opposite Jz col- 
liding, giving two paraexcitons with momenta K' and P'. 
Fermi's golden rule gives for the rate 



K,P,K',P' 

x5{Ek' + Ep, -Ek-Ep- 2AE)Sk+p,k'+p', (1) 

where M is the matrix element for this process, is the 
distribution function of the i species (ortho or para exci- 
tons), having a dispersion relation i^K = h^K^/2m, with 
m being the total exciton mass. In this crude calculation 
we consider a cold (K, P w 0) orthoexciton gas, which 
allows us to write that 

ro,p«-7v2 \M\\i + f^,)ii + f^,) 
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(2) 



where Ni is the total number of excitons of species i. The 
energy conservation condition in the above equation im- 
plies that K' and P' are of order {m/S.E /h?)^/"^ . Since 
for these wavevectors the occupation number is much less 
than 1, we can ignore the enhancement factors 1 + / 
above. In addition, we argue below that the typical mo- 
mentum exchange that enters the matrix element M is 
of order a^" - {mEb/h^y/'^ > {mAE/h^y^'^ , where as 
is the exciton Bohr radius, and Ej, is the exciton binding 
energy. Since Ei, 3> AE, M does not vary substantially 
in the sum of Eq. (0) and can be taken outside it. 
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The last sum is simply the density of states calculated for 
an energy AE. The interaction that enters the matrix 
element M is some screened Coulomb potential y(g, w). 
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where f2 is the volume of the crystal, and e(q, a;) is the 
dielectric function . The exciton wavefunction $k (re — 
r^) of an exciton carrying momentum K can be written 
as 



^'K(re 
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„j:K-(r,+rh)/2 
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where and are the electron and hole coordi- 
nates, and 0q 8(7ra^)^/^/[l -|- (qas)^]^ is the Fourier 
transform of the ground state hydrogenic wavefunction, 
$ = e-'/"'' /{Tra%), which we assume to be the relative 
electron-hole wavefunction. We have assumed that the 
two colliding orthoexcitons have K = P = 0; denoting 
the momenta of the electron and the hole in each pair as 
k, — k, and p, — p [see Fig. 1], after the two excitons have 
exchanged their electrons or holes, and some momentum 
q, conservation of energy and momentum requires that 
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= AE. 



(6) 
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FIG. 1. The process of two orthoexcitons with opposite 
J2 colliding, exchanging their holes and some momentum q, 
giving two paraexcitons in the final state. Time progresses 
from left to right. The straight lines with arrows pointing to 
the right (left) denote electrons (holes). The wiggly line de- 
notes the Coulomb interaction. The arrows pointing upwards 
(downwards) denote Jz = 1/2 {Jz = —1/2) electrons or holes. 



Since the exciton wavefunction has a momentum spread 



of order , we see that p ~ fc ~ , or in other words, 
?i^p^/2m ~ Eb. Because of the condition /S.E <C i?b, 
we conclude that q is also of order a^^ Further- 
more, it is a rather good approximation to assume that 
e(q, w) ~ Co [^i the low-frequency dielectric constant of 
CU2O, which is approximately equal to 7.5, and thus 
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After some trivial manipulations we can write for the rate 
of the orthoexciton conversion process Fo.p/ri = no/To p, 
where Ui is the density of species «, and 



(8) 



Let us now examine the reverse process of two paraex- 
citons colliding, giving two orthoexcitons. The energy 
scale which is very important for this process is the en- 
ergy splitting AE between the orthoexcitons and paraex- 
citons in CU2O, which is equal to 12 meV |^ at the zone 
center, corresponding to a temperature of approximately 
150 K. We can write for the decay rate Fp o/f2 — Up/Tp o, 
with 



To,p {np/no)e 



where T is the temperature of the exciton gas. Thus, for 
temperatures which are much lower than AE jkB-, the in- 
terconversion mechanism converts the orthoexcitons into 
paraexcitons, but not the reverse, since the para-ortho 
process is exponentially suppressed [e-^^/'^s'^ 0]. On 
the other hand, for temperatures which are > AE/ks: 
the two rates Fo,p and Fp ^ can be comparable, depend- 
ing on the relative densities of ortho and paraexcitons, 
and thus the net rate can be very low. 

To calculate the actual value of the rate given by 
Eq. (|§), we use the following numbers for excitons in 



CU2O (which have very low uncertainty): the binding 
energy is 153 meV, the energy splitting AE at the zone 
center is 12 meV, and finally the Bohr radius as is 5.3 
A, as given by a variational calculation presented in Ref. 
P|. With these numbers we get 
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) ns 
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where the notation n(10 cm^'') means that the density 
is to be measured in units of lO^^cm"'^. 

We can now compare the experimental decay rate of 
orthoexcitons that was measured recently by O'Hara et 
al. [01 with the above theoretical number on the one hand, 
and with the theoretical number for the Auger process on 
the other hand. For the low lattice temperature of ~ 2 K 
in this experiment, the average kinetic energy per particle 
of the gas of orthoexcitons is expected to be much less 
than AE, and thus we assume for the net interconversion 
rate F^ p Tp q ~ Fq p. The "two-body decay constant" , 
A K, 10~^^ cm'^/ns, that is extracted from the data given 
in Eq. (5) of the first paper in Refs. § has the same order 
of magnitude as the rate given by Eq. (|lO|) . 

We mentioned before that the Auger process is ex- 
pected theoretically to have a rather low decay rate, in 
view of the very long radiative lifetime of orthoexcitons 
which was reported recently in Refs. ||l[. More specifi- 
cally, in the theoretical study of the Auger process de- 
scribed in Ref. |11|, a detailed analysis of this mecha- 
nism implied that the phonon-assisted Auger decay pro- 
cess was the dominant one. However, to get the rate, the 
authors used the orthoexciton radiative lifetime of 25 ns 
at a temperature of 10 K |l^ that was measured in the 
past. In Refs. ||l| the same quantity was measured to 
be approximately 10 yits, and since the radiative lifetime 
can be limited by imperfections or any other factor, the 
real radiative orthoexciton lifetime is at least 10 /is, or 
even longer. This implies that the phonon-assisted Auger 
decay rate, based on the theoretical study of Ref. Q is 
negligible. These arguments provide strong evidence that 
the ortho to para exciton interconversion mechanism that 
we study here is really the dominant process. 

We turn now to the contradiction between the two 
methods which have been used for determining the or- 
thoexciton density. It is important to get first an es- 
timate for the interparticle clastic scattering rates, in 
order to compare them with the rates of interconver- 



(9) sion by spin-exchange. Typical scattering rates r 



for elastic collisions between excitons are expected to 
be on the order of = ncrwth, where u is the scat- 
tering cross section, and Wth is the thermal velocity, 
"^th = {^ksT / 1:111)^/ . At the low temperatures of inter- 
est one can assume hard-sphere scattering between the 
excitons. If a is the scattering length, then for identical 
bosons a = Stto^. Recently the scattering length for ex- 
citons in CU2O has been calculated to be on the order of 
2aB [0 with use of Monte-Carlo simulations. Therefore 
we get 



0.1n(10^^cm"3) VT(K): 



(11) 
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assuminff that the total exciton mass is equal to 3 electron 
masses 1^. Here the density is measured in units of 10^® 
cm~^, and the temperature in units of degrees Kelvin. 

After the laser pulse starts to decrease, the paraexci- 
tons become the dominant component of the gas because 
of the interconversion mechanism. For typical paraexci- 
ton densities of 10^^ cm~^ and thermal velocities of order 
50 K, one sees from Eq. ( ITll ) that the typical scattering 
times are of order 100 ps, and the paraexcitons should 
be able to establish thermal equilibrium. The paraexci- 
tons should also have a well-defined chemical potential, 
since after the initial times that the orthoexcitons con- 
vert to paraexcitons due to the interconversion process, 
their number does not vary significantly with time. By 
contrast, the orthoexciton-orthoexciton elastic scattering 
processes become less and less frequent because of their 
decreasing density, even if the orthoexciton-paraexciton 
collisions can bring them to thermal equilibrium; we 
claim, however, that chemical equilibrium has not been 
established in the orthoexciton gas, since the orthoexci- 
tons have a relatively fast way of converting into paraex- 
citons and their number is not conserved. One can spec- 
ulate that under such circumstances the orthoexciton gas 
could have a chemical potential which is rather low, but 
this is a non-equilibrium problem, and it requires a de- 
tailed study. One, for example, could use the Boltzmann 
equation to describe all the important processes which 
take place, and derive the distribution function of the 
orthoexcitons as a function of time. 

There are some remarks one can make concerning this 
model we are proposing. Firstly, the estimate we made 
for the interconversion rate, Eqs. (||) and (||), does not 
assume thermal equilibrium, which is rather important in 
our problem. Secondly the Bose-Einstein condensation of 
paraexcitons does not seem hard, since as we explained 
there is no efficient mechanism which would destroy them 
on the timescales of interest, and their expansion could be 
the only factor against them in order to Bose condense. 
However the expansion is not expected to be dramatic, 
and it can also be reduced by applying some stress to the 
crystal, thus effectively trapping the excitons. 

Finally, it has been argued is Ref. , that the orthoex- 
citons - provided they are not far away from equilibrium 
- have to move along lines of constant entropy, along 
which Uo oc T^/'^. To derive this result, the authors 
assumed that there is a competition between acoustic- 
phonon cooling and Auger heating. Since the phonon 
cooling rate is oc —T^l"^ for low lattice temperatures, and 
the Auger heating rate is esc rto, no (x r3/2. Remark- 
ably, if indeed the Auger process is negligible, and the 
interconversion process is the dominant mechanism, the 
heating rate due to this effect is equal to Ai?/ro^p, and 
thus still proportional to tIq. Since this argument does 
not depend on the quantum degeneracy of the g as, tIq 
should still be proportional to T^l^: we conclude again 
that the orthoexcitons are expected to move parallel to 
the phase boundary, along adiabats, in contrast to the 
paraexcitons which most probably Bose condense. More 



experimental and theoretical work are required to verify 
these predictions. 
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